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ABSTRACT 

A liquid mixture viscosity as a function of composition is extremely complex. Theoretical considerations have 
offered little help in explaining these complexities. Attempts made so far to derive a generalized expression for 
viscosities of all mixtures resulted in equations with many undetermined constants, and no method allows a 
reliable prior prediction of these constants. These methods, therefore, are purely descriptive. In this paper, we 
made an attempt to analyse the viscosity of liquid mixtures employing Schrddinger equation. We came up with a 
model >// . ( r ) y/ . ( r ) = ^ y/ n . ( r ) y/ n . ( r ) exp (E r . /k B T) that should account for the viscosity coefficient 

of a liquid mixture when the knowledge of the potential U„(r) + XUj.„(r) would be available for the molecules of 
the liquid mixture. 
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I. INTRODUCTION 

Some pairs of liquids are miscible in all proportions and complete series of mixed crystals are formed 
by certain pairs of crystalline solids. Other pairs of liquids and solid substances are practically immiscible. A 
third group is formed by pairs which are miscible in limited proportions. Miscibility depends upon the 
temperature, and in general is favoured by the chemical similarity of the substances concerned. Water does not 
mix with hydrocarbons, is miscible with methyl and ethyl alcohols, partially miscible with butyl and amyl 
alcohols, but practically immiscible with higher members of the alcohol series. When two partially miscible 
liquid substances A and B are brought together, then, in general, two layers are formed which contain different 
proportions of A and B. The composition of these layers so-called conjugate pairs, (Wohlfarth, 2009) depends 
on the temperature. As a rule, the difference between the conjugates diminishes with rise of temperature and in 
some cases may ultimately disappear. The temperature at which this occurs is the critical solution temperature 
and above this the two substances are miscible in all proportions. Increased miscibility is sometimes produced 
by a fall of temperature and a few pairs of substances are known for which both upper and lower critical 
temperatures can be observed. Nicotine and water, for example, are completely miscible below 61° C and above 
210° C; between these limits they are only partially miscible. 

Since the properties of pure liquids are largely dependent on the attractive forces between the molecules, it is 
extremely unlikely that the relations connecting the properties of liquid mixtures with those of their components can be of 
the same simple nature as those which characterize the corresponding gaseous systems. In some cases the mixing of two 
liquids is accompanied by very pronounced changes in volume, by a rise or fall of temperature and by other effects. For 
other liquid pairs, such effects are so small as to be scarcely measurable. The relations between the properties of liquid 
mixtures and those of their components may be conveniently shown by plotting the measured property against the 
composition of the mixture. The property-composition curves obtained in this way sometimes deviate but little from the 
straight line which corresponds with the simple mixture rule, but for other mixtures large deviations are found, and 
frequently the curves show a maximum or a minimum, (Hulya, 2000). These deviations result to the excess values which are 
the differences between the experimental and the theoretical values. 

II. METHOD 

For a temperature range from the freezing point to somewhere around the normal boiling temperature, it is often a 
good approximation to assume that the logarithm of the viscosity of liquids is linear in the reciprocal absolute temperature. 
In general, pure liquid viscosities at high reduced temperatures are usually correlated with some variation of the law of 
corresponding states, such as the model by Sastri (1992). At low temperatures, most methods are empirical and involve a 
group contribution approach. Current liquid mixture correlations are essentially mixing rules relating pure component 
viscosities to composition. Little theory has been shown to be applicable to estimating liquid viscosities within a reasonable 
accuracy. 
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Almost all methods to estimate or correlate liquid mixture viscosities assume that values of the pure 
component viscosities are available. Thus the methods are, in reality, interpolative. Nevertheless, there is no 
agreement on the best way to carry out the interpolation. Irving (1977) surveyed more than 50 equations for 
binary liquid viscosities and classified them by type. He pointed out that only very few do not have some 
adjustable constant that must be determined by experimental data and the few that do not require such a 
parameter are applicable only to systems of similar components with comparable viscosities. He recommended 
the one-constant Grunberg-Nissan equation as being widely applicable with reasonable accuracy except for 
aqueous solutions (Grunberg and Nissan, 1949 

The predictive methods for viscosity of liquid mixtures include semi -theoretical and empirical models. 
Most semi-theoretical models (Shuai et al, 1984; Barrufet and Setiadarma, 2003; Al-Besharah et al, 1987; Liu et 
al, 1999 and Li, 1992) for petroleum fractions, which have a theoretical framework but parameters determined 
from experimental data, are based on either the corresponding-states approach or the modified Chapman-Enskog 
theory. A wide variety of liquid-mixture viscosity prediction formulas (Irving, 1977) use the simple-mixing-rule 
equation based on calculations of the weighted average of the component viscosities. In order to also contribute 
by different approaches to get reasonable theoretical model for the analysis of viscosity coefficient of liquid 
mixtures, we decided to do our approach from the quantum point of view. 

III. RESULTS AND DISCUSSIONS 

The viscosity coefficient, u of a liquid, (De and Dikko, 2012) is given by 

p = B exp (Ea/k B T) (1) 
For a multi component liquid mixtures 

P = PiMi + P 2 M 2 + + P„P„ = X P nV„ (2) 

where p a is the viscosity coefficient of the liquid component n, whose mole fraction in the mixture is P„, and P„ 
is related to Vi,V 2 V„ by the equation 

(VoJm ) 

P = v - HhI n - (3) 

iy r p i ju , + v 2 p 2 /m 2 + +v„p„ /m n ) 

where p n and M n are the densities and molecular weights of the components respectively. The experimentally 
determined, u cx will differ from u of equation (2) by 

p a = p + p E (4) u E is the 

excess viscosity which is a correction to equation (2) arising out of interaction that include clusters formation, 
caging, etc, (Nagasawa et al, 2003) among the molecules of different components of the liquid mixture. 
Assuming that the interaction between molecules are of a pair wise nature accomplished by central forces 
having a potential £/(r) and a radial distribution G(r). These two important parameters determine the equation of 
state of the liquid. If the liquid is assumed to consist of spherically symmetric molecules the equation of state is 
given by the following equation: 

—G(r)r 2 dr (5) 

[ 3VK J y° dr 

where P is the pressure, Temperature, K B is the Boltzmann constant and V is the average volume per particle of 
liquid. 

In a pure liquid, the distribution function G n (r) may be written as 

G ii (r) = Y J P„ i ¥ ffr)¥ „(r)' ( 6 ) 
!F n is the quantum mechanical wave function of a molecule of the liquid type n at position r and the Schrodinger 
equation for the molecule of nth component liquid may be expressed as 

(-h 2 v 2 \ 

\~i +f/ »<» k»<» = E nV„(r) (7) 

P n is expressed in average number of molecules of the nth type per unit volume of the liquid. In a multi- 
component liquid mixture, 

G(r) = £ G n (r) and [/(O^^.W + ^f/^W (8) 

Where 
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The overall Schrodinger equation for the mixture of liquids may be written as 

r n 2 v 2 i 

+U (r) pFO) = 



(10) 



(9) 



£/(r) is the potential at a point r arising out of interaction between molecules in the liquid mixture. The exponent 

" of p n in the expression for ^(r) in equation (9) is based on the assumption that the probability of finding a 
molecule of the nth type liquid in a small volume AV of the liquid mixture about the position r should be 
proportional to p n V and not to any other power of p n V. The average energy per particle can then be written as: 



G„(r) depends on P„ of equation (2). Ideally it is expected to have also an implicit dependence on temperature also, at least 
through the term P. E is the activation energy embedded in the second term of the above equation (11) for the average 
energy. Once E a can be evaluated, the [i of the resulting liquid mixture and hence (j E can be evaluated from equations (1) and 
(11). This correction of p E is obviously a function of P n , true nature of which is not yet theoretically known but expected to 
be obtained from above equations. 

The detailed explanation for this very interesting observation from the molecular physics point of view is left for 
future work including the theoretical computation of p E from above model. Once p E can be understood theoretically then 
equations (2) - (4) can form the basis of a new emerging technology of analyzing the concentration of the liquid components 
of any liquid mixture which does not contain dissolved solid materials. 

From equation (11) the temperature variation of activation energy E a can be accounted for in the first order by 
considering the thermal expansion of the volume and consequently the thermal variation of G(r) through the term P. It may 
be mentioned that the thermal expansion coefficients (Cutnell et at, 1995) for methyl alcohol, ethyl alcohol and water 
respectively: 1240, 1200 and 207 times 10" 6 . From above theory it is expected that -6E a /5T for ethyl alcohol would be about 
six times than for alcohol-water mixtures. This latter theoretical conjecture yet remains to be validated by experiments. One 
thing that is clear from the above model is that it can clearly account for the decrease of \i of liquid with temperature. 
Another dependence on temperature may also come from the Boltzmann population of the different eigen energy levels E of 
equation (6) in which case ^rj^r)* in equation (6) may be replaced as 



When the knowledge of U„(r) + SUi_ n (r) would be available for the molecules of the liquid mixture, the above 
model together with equation (1) should account for the viscosity coefficient of a liquid mixture. 

IV. CONCLUSION 

The model developed in this paper reveals that with further studies and the knowledge of the potential 
for the molecules of the liquid mixture, a better theoretical model can be developed for the analysis of multi- 
component liquid mixtures. 
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